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Abstract 

We consider an ideal mesoscopic ribbon in which a steady azimuthal current 

is generated. We show that the closed interacting electron system in the 

presence of the current is described by a density matrix which is that of an 

equilibrium system without current but with a constrained Hamiltonian. 
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Two recent papers present new approaches to nonequilibrium steady-state systems. The 
important common point made in both papers is that steady-state mesoscopic systems can 
be described by a density matrix which has the form of an equilibrium density matrix, but 
with a constrained Hamihonian. In the first paper, Hershfield |I| demonstrated that the 
density matrix p of a steady-state nonequihbrium quantum system has the general form 

p = e-^(^"^), (1) 

where H is the Hamiltonian, P is the inverse temperature (we will use units in which fc^ = 
h = e = nie = 1), and Y is an operator which depends on how the system is driven 
out of equilibrium. This operator is defined implicitly in terms of a non-trivial infinite 
set of differential equations. In the other paper (Ref. 0], hereafter referred to as HJ), we 
formulated an approach to steady-state mesoscopic transport based on the maximum entropy 
principle of nonequilibrium statistical mechanics . Using the maximum entropy principle, 
we derived the following density matrix for a multi-terminal steady-state mesoscopic system: 



where is the current incoming from terminal a and the Lagrangian multipliers are 
adjusted to give the applied source-drain current and zero current at all other terminals. 
This density matrix, like that in Eq. (|l]) , has the form of an equilibrium density matrix of a 
constrained Hamiltonian. It can also be argued on general grounds using Galilean invariance 
that distributions of the form obtained from Eq. (Q) can be expected in an infinitely long 
ideal mesoscopic wire 

One important point to be noted is that the single-particle distributions obtained in HJ 
differ, even when linearized with respect to the currents, from the local-equilibrium 
distributions typically used in the Landauer-Biittiker formalism of mesoscopic transport. 
In the linear-response regime, one would not expect that measurements of, for example, 
conductivity directly probe the electron distributions, since such quantities can in the linear- 
response approximation be expressed as the trace over the equilibrium density matrix and 




(2) 
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two-point correlation functions. However, experiments conducted beyond the linear-response 
regime should be sensitive to the nonequilibrium steady-state electron distributions. This 
should be the case for precision measurements of the integer quantum Hall effect, where the 
measured Hall voltage exceeds IQTiojc/e, with ojc the cyclotron frequency. We demonstrated 
in HJ that a transport theory based on Eq. (j^) can explain the quantization of the integer 
quantum Hall effect at such large, but experimentally typical, currents, while the Landauer- 
Biittiker formalism fails to do so. 

In deriving the density matrix Eq. (^) using the maximum entropy principle, it was 
necessary to make assumptions about which are the relevant observables, the expectation 
values of which are taken to be known. While this is a standard procedure of the maximum 
entropy approach to statistical mechanics, it is also the source of some controversy about 
this approach, since no unambiguous procedure for choosing the observables exists 0. It is 
therefore desirable to inquire whether such steady-state distributions in mesoscopic systems 
can be obtained by other means. The purpose of the present paper is to present one such 
example. We will show by an explicit calculation that the exact density matrix in the 
presence of a steady current is precisely that of HJ for a simple specific case: a closed system 
of interacting (spinless) electrons in an ideal, mesoscopic, two-dimensional ribbon. Here we 
obtain this result by considerations of equilibrium thermodynamics in a rotating reference 
frame, plus adiabatic switching-on, without appeal to the maximum entropy principle. By 
an ideal system we mean that there are no elastic or inelastic scattering processes other than 
those resulting from electron-electron interactions. A current is generated in the electron 
system by adiabatically threading the ribbon with an integral number of flux quanta. In 
this case, there is then a one-to-one correspondence between the many-body eigenstates 
in the presence and in the absence of the current. This allows us to demonstrate that 
the density matrix in the presence of the current can be directly related to the density 
matrix in the absence of current, and has the form of an equilibrium density matrix with 
a constrained Hamiltonian. While we draw no conclusions about open systems here, we 
provide an explicit example of a system which by an independent calculation proves to be 
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described by the density matrix obtained by HJ. 

The system considered here resembles those considered in investigations of persistent 
currents [§] in mesoscopic rings. The important difference is that in those investigations 
the current as a function of an apphed AC magnetic flux is typically calculated for a non- 
interacting system in the presence of elastic scatterers. (The role of interactions in the 
presence of disorder is complicated and unclear. |^) Great care has to be taken to ensemble- 
average correctly, and to correctly account for the magnetic field actually penetrating the 
ring itself in an experiment. Also, in the presence of scatterers which break the rotational 
invariance there is no adiabatic curve-crossing as the magnetic flux is increased adiabatically. 
Here we are considering a closed, impurity-free interacting electron system, and the flux is 
used only as a device [|1^ to generate the electric field and the resulting current. Our 



ultimate goal is to obtain the exact density matrix in the presence of a steady current. 

We take the system to be a two-dimensional interacting electron gas confined to a cylin- 
drical ribbon of radius R. Positions on the ribbon are described by cylindrical coordinates 
(r, 6, z), with r = R. A uniform azimuthal electric field E6 is generated by piercing the bore 
of the cylinder with adiabatically increasing magnetic flux. The electric field is described by 
a time-dependent vector potential 

A{r,e;t) = ^f{t)e, (3) 

where $/ is the magnetic flux piercing the cylinder's bore. The monotonically non- 
decreasing function f{t) describes the adiabatic turning-on of the vector potential, with 
f{t — oo) = and f{t oo) = 1. We will assume that $ = p^o with p an integer and 
$0 = 27rc the flux quantum, so that the cylinder contains an integral number of flux quanta 
as t oo. From 

the electric field then vanishes for t ±oo but is finite during the time that / is changing. 
This finite electric field sets up an azimuthal current which persists as the electric field 
vanishes, since the system is dissipationless. 



With the vector potential given by Eq. (j^), the first-quantized many-body Hamiltonian 



IS 



2m* 



(5) 



where the sums are over the particles of the system, Vc{z) is a confining potential, Vij = 
|rj — Tjl, and V{r) is the electron-electron interaction. Note that j = —id/dOi is the 
canonical angular momentum operator, which is independent of A. The eigenstates of H 
are given by A^-particle wavefunctions \E' = ^{ziOi, Z2O2, • • • , znOn)- To proceed, we second- 
quantize the Hamiltonian Eq. (|^) using the field operators 



1p{z,9) = '^Cmni^mniz, 
ran 

where the single-particle wavefunctions ipmniz^O) = 
probability per unit area) satisfy 
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/2-kR 



2m 



+ V,{z) 



(6) 



z) (normalized to give 
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Here n is a subband index and m an angular momentum index. The second-quantized 
Hamiltonian is 
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[pf{t)f + 2mpf{t) 
2m*R^ 



mn mn 



mm' k 



(8) 



"1"2"3"4 

with Vfc;ni,n2,n3,7i4 the matrix elements of the particle interaction. The total angular momen- 
tum operator = J2i Lz,i becomes 



(9) 



The azimuthal current density operator is jg{z,9) = ^[ip'^vg'tp + {vgtp)^tp], where vg 
2m*R {.Lz + pf) is the speed in the azimuthal direction. This can be written 
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2'Km*B? 
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m' + m 



Pf{t) 



(10) 



We also introduce an azimuthal current operator / by integrating jg{z,6) over z and aver- 
aging over 6, which yields 

i 1 



5][m + p/(t)]4„c^„. (11) 



It is straightforward to show that H, /, L^, and the number operator N = J2mn '^nn'^mn 
commute, so we choose a basis in which all are diagonal. In fact, 

i 1 



;{L,+pfN), (12) 



so a basis diagonal in and N is automatically diagonal in /, and the eigenvalues I and 
M of I and are simply related. 

Let us consider next the adiabatic evolution of some iV-electron wavefunction \E'. Suppose 
that is an eigenstate of Hq = H{f = 0) with energy E and total angular momentum M. 
From Eq. (^, the Hamiltonian in the presence of the flux tube can be written 

As long as the field is switched on adiabatically {i.e., df /dt is much smaller than electron- 
electron scattering rates), \Ef remains an eigenstate of H when / 7^ 0. Although its (canoni- 
cal) total angular momentum is unchanged, its energy and current change with /. Since p is 
an integer, the vector potential as / ^ 1 can be removed by a unitary gauge transformation: 

H = e'P^^^^He-'P^^^^ = Hq. (14) 

The transformed state has total angular momentum M + pN, and is an eigenstate of Hq. 
Thus, the spectrum of the Hamiltonian for t ^ 00 is identical to the spectrum for t —>■ —00. 
More, the energy and current when / = 1 of a state with angular momentum M are precisely 
those of an eigenstate of Hq with angular momentum M + pN. It is clear from this and 
Eq. (^) that as / grows from to 1 the subset of the spectrum consisting of eigenvalues 
of states with angular momentum M evolves to the subset of eigenvalues of eigenstates of 



Hq with angular momentum M +pN. This occurs with no level crossings between states in 
this subspace (although crossings occur between states with different M). 

In particular, suppose we start at t ^ — oo with the system in the interacting ground 
state (with M = 0) and adiabatically turn on the vector potential. Then at t ^ oo the 
system's energy and current will be that of the lowest-energy eigenstate of Hq with angular 
momentum pN. Thus, the final energy and current are obtained by finding the lowest-energy 
eigenstate of Hq in the subspace with (Lz) = pN. In other words, the final state when a 
current has been turned on can be found by extremizing Hq subject to the constraint that 
= pN. Using Eq. (0), we can instead constrain the current to be (Jq) = pN/2nm*R'^ 
(where Iq is the current operator with / = 0). In practice, a convenient way to satisfy the 
constraint is to introduce a Langrangian multiplier, and to look for stationary states of 

Ho - e/o (15) 

where ^ is chosen so that Iq has the required eigenvalue. 

We next turn to finite temperatures, and consider a system initially in contact with a 
particle reservoir, which maintains the average number of particles at A^, and in thermal 
contact with a heat reservoir at temperature T. Initially (t —00) the density matrix is 

and the system has zero average current, (/o)t = Tr{lQpQ)/Tr{pQ) = 0. We then insulate 
the system thermally from the heat reservoir and adiabatically turn on the vector potential 
while keeping the average number fixed. There are two ways to understand what then 
happens. First, notice that when the system is isolated thermally each subspace of states 
with a given total angular momentum M becomes a separate subsystem. Although electron- 
electron interactions could cause re-equilibration within a subset, these interactions conserve 
total angular momentum and so cannot, e.g., transfer energy from one angular momentum 
subspace to another. Moreover, as explained above, the energy of each state in a given 
subset changes by exactly the same amount when / varies (see Eq. (|T3|) ). That is, the 
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spacings in energy between the states in each subspace remain constant, and consequently 
the occupancies of these states will not change. Then the occupancy of a state at t — > oo is 
given precisely by its initial occupancy at t ^ — oo. Consider an A^-electron state with initial 
energy E and angular momentum M. Initially the occupancy of this state is e~^(^~^o^). At 
t -—' CO, this energy and current of this state become 

= ^ + + 

The occupancy of this state is still e~^^^~^"^\ That is, the occupancy is determined not by 
its energy E', but by E, which can be related to E' and /' using the above. Since the initial 
and final sets of states are identical (for p an integer), the density matrix can be written in 
terms of the / = operators Hq and Jq. As t oo, the density matrix evolves to 

p = exp{-a(ff„--Ji^(L, + ipiv)-,„iv)} 
= exp{-/j(//„-a„-fiA')}, (18) 

where the final equality defines ^ = 2ttp and /i = /iq + p'^/2m*R'^. 

A second way to understand this result is to consider equilibrium in a rotating reference 
frame. As the vector potential is turned on, the system evolves adiabatically as t ^ oo 
to a state with average current {I)t = I- After the electric field has returned to zero (as 
t — > cx)), we transform to a coordinate frame rotating with an angular velocity Vt relative to 
the lab frame, where VL is chosen so that the azimuthal current is zero in the rotating frame. 
The Hamiltonian H in the rotating frame and the Hamiltonian H in the original frame are 
related by [|n|]: 

H = H- ni,, (19) 

where 

n = (20) 
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In the rotating frame, the system consists of an isolated interacting electron gas at zero 
net current. In this frame, one should expect the electron-electron interactions to equilibrate 
the system. (This equilibration is at the same temperature T that the system was in at 
t —oo in the lab frame, since the center-of-mass motion of a system does not change 
the temperature |T2[-) Thus, in the rotating frame, the system has an equilibrium density 



matrix 



p = e 



-/3(^-A7v)_ ^21) 



We express this density operator in terms of the Hamiltonian Hq of the stationary frame at 
zero flux by first using Eq. (|19D and then gauging away the vector potential as above. With 
the operators in their first-quantized form, this gauge transformation gives 

where n = fl + p'^ / m* R"^ . Thus the system is described by the density matrix 

^ = e-/3(i?o-?/o-MJv)^ (23) 

where we have used Eq. (|I2]) (with Jq = /(/ = 0)), and as before ^ = 2np. In practice one 
can regard fi and C, as parameters which are adjusted to give the required thermal averages 
of particle number and current. Thus the density matrix is that of an equilibrium systems 
described by a constrained Hamiltonian, and is identical to that (Eq. (^) obtained in HJ 
by maximizing the entropy subject to constraints on internal energy, particle number and 
total current. 

This work was supported by the National Science Foundation through grant DMR93- 
01433. 
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